A statistical model analysis of yields and fluctuations in 200 GeV Au-Au collisions 
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We show that the simultaneous measurement of yields and fluctuations is capable of falsifying 
and constraining the statistical hadronization model. We show how such a measurement can test 
for chemical non-equilibrium, and distinguish between a high temperature chemically equilibrated 
freeze-out from a supercooled freeze-out with an over-saturated phase space. We perform a fit, and 
show that both yields and fluctuations measured at RHIC 200 GeV can be accounted for within the 
second scenario, with both the light and strange quark phase space saturated significantly above 
detailed balance. We point to the simultaneous fit of the K/tt fluctuation and the K* /K~ ratio as 
evidence that the effect of hadronic re-interactions after freeze-out is small. 

PACS numbers: 25.75.-q,24.60.-k,24.10.Pa 



The statistical hadronization model (SHM) 0, 0, Q 
has been extensively applied to the study of soft particle 
production in hadronic systems. When it includes the 
full spectrum of hadronic resonances Q, the SHM can 
describe quantitatively the abundances of all light and 
strange hadrons produced in heavy ion collisions. 

The ability of the SHM to describe not just averages, 
but event-by-event multiplicity fluctuations has however 
not been widely investigated, and its applicability is cur- 
rently a matter of controversy. Event-by-event particle 
fluctuations have elicited theoretical |E S 13 an d ex- 
perimental interest M, yUj , both as a consistency check 
for existing models pa, ||J , as a test of hadron gas equili- 
bration and the presence of resonances [j| , and as a way 
to search for new physics 0. 

This study investigates weather the SHM can explain 
both yields and fluctuations with the same parameters. 
We illustrate the use of both yields and fluctuations as 
a probe capable of differentiating between freeze-out sce- 
narios. We then perform a fit to preliminary 200 GeV 
RHIC experimental data, and discuss the results in the 
context of chemical (non) equilibrium within the RHIC 
system, the freeze-out temperature and the presence (or 
absence) of hadronic reinteraction after freeze-out. 

The statistical hadronization model assumes that par- 
ticles are produced according to a probability determined 
by their phase space density. The first and second cumu- 
lants of this probability distribution give, respectively, 
the average value (over all events) of the desired observ- 
able, and its event-by-event fluctuation. 

Conserved quantities can be treated in several ways, 
appropriate to different experimental situations: If the 
totality of the system is observed, than conserved quan- 
tities can not fluctuate. If a small fraction equilibrated 
with the rest of the system ( the "bath") is observed, 



than conserved quantities will fluctuate, event-by-event. 
Rigorous conservation is known as the Microcanonical en- 
semble, while allowing energy and other conserved quan- 
tities to fluctuate between the system and the bath leads, 
respectively, to the Canonical and Grand Canonical (GC) 
ensembles. As shown recently all fluctuations are 
ensemble-specific even in the thermodynamic limit. 

In this work, we use the GC ensemble, implemented 
in a publicly-available software package to calcu- 
late fluctuations and yields. We motivate this choice by 
the fact that RHIC experiments observe the mid-rapidity 
slice of the system, comprising roughly 1/8 of the total 
multiplicity, an appropriate fraction for a GC prescrip- 
tion. Boost invariance, a good symmetry around mid- 
rapidity, links this rapidity slice with a sector in configu- 
ration space. If the system observed at RHIC is a nearly 
ideal fluid, the matter created in this space should be 
in equilibrium, grand-canonically, with the unobserved 
regions. If freeze-out temperature throughout observed 
space is approximately constant, the GC ensemble should 
be able to describe both yields and fluctuations [Hi . Il4| . 

The final state yield of particle i is computed by adding 
the direct yield and all resonance decay feed-downs. 
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Here (Ni) is the ensemble average of the i-th particle 
(resonance) multiplicity at chemical freeze-out and Bj— >j 
is the branching ratio of the j — » i decay. The fluctuation 
(AX 2 ) = (X 2 ) - (X) 2 also has two parts 

(AN 2 ^) = B j ^ i (l-B j ^ i )(N j )+B^ i (AN^) (2) 

assuming no primordial correlation between i and j. 

The fluctuation of a ratio N\ /N% (a useful experimental 
observable, since volume fluctuations cancel out event- 
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by-event) can be computed from the fluctuation of the 
denominator and the numerator Q [a\ = (AX 2 ) / (X)): 
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Note the appearance of a negative correlation term Di 2 , 
which arises due to the correlation between N% and N 2 
stemming from a common resonance feed-down (A — > pit 
will be a source of correlation between N p and iV w ). 

The average yields and fluctuations of hadrons them- 
selves can be calculated by a textbook method given a 
grand canonical description of the phase space; For a 
hadron with an energy E p = \/p 2 + m 2 , the average 
number, and event by event fluctuation will be given by 
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where the upper sign is for Fcrmions and the lower sign 
is for Bosons. Here A; is the fugacity of particle i, see 
Eq. I jlOj l. These formulae can be efficiently calculated 
through a Bessel function expansion 0, 0] . 



Note that the series start diverging at Xi 



and that for fluctuations the higher order terms in the 
expansion are always more significant than for yields: for 

A mJT =1 _ e 
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(9) 



The parameter A; corresponds to the particle fugacity, 
related to the chemical potential by Xi — e Mi / T . Pro- 
vided the law of mass action holds, it should be given by 
the product of charge fugacities (flavor, isospin etc.). ft 
is then convenient to parametrize it in terms of equilib- 
rium fugacities A° q and phase space occupancies 7^. For 
a hadron with q(q) light quarks, s(s) strange quarks and 
isospin I3 the fugacity is then 
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If the system is in chemical equilibrium then detailed bal- 
ance requires that j q — "f s — 1. In an expanding system, 
however, the condition of chemical equilibrium no longer 
holds. Kinetically, this occurs because collective expan- 
sion and cooling will make it impossible for endothermic 
and exothermic reactions, or for creation and destruction 
reactions of a rare particle, to be balanced. Provided the 
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FIG. 1: Fit of preliminary 200 GeV data, including the 
A' ± /7r ± fluctuations and the K* (892) and A* (1520) resonance 



system remains in local thermal equilibrium, Xi can still 
be used as a Lagrange multiplier for the particle density, 
and the first and second cumulants can be calculated 
from the partition function the usual way [TH Il6| . How- 
ever, in this case in general j q 7^ 1,7s 7^ 1. 

This picture becomes particularly appropriate if the 
expanding system undergoes a fast phase transition from 
a QGP to a hadron gas. In this case, chemical non- 
equilibrium and super-cooling 0] can arise due 
to entropy conservation: By dropping the hadroniza- 
tion temperature to ~ lAOMeV and oversaturating the 
hadronic phase space above equilibrium (7^ ~ 1.5, 7 S ~ 
2), it is possible to match the entropy of a hadron gas 
with that of a system of nearly massless partons • 

Fits to exp erimental data at both SPS 0] and RHIC 
energies |20| indeed support these values of 7 9 . s when 
these parameters are fitted. Moreover, best fit 7 9)S > 1 
arises for a critical ener gy |2l| (correspondin g to the en- 
ergy of the K/n "horn" 22] ) and system size [23J, just as 
expected from the interpretation of j q as a manifestation 
of a phase transition. However, the fits performed in |2*l) 
have not been able to rule out equilibrium models (at 
SPS and RHIC the difference in statistical significance 
between equilibrium and non-equilibrium is ~ 20%). 

In the absence of contrary evidence, equilibrium mod- 
els are usually preferred for their smaller number of free 
parameters [24. |25| . Equilibrium freeze-out temperature 
varies between fits, ranging from 155 [2lJ to 177 [25| MeV. 

Both scenarios are physically reasonable, can describe 
the data, and would be instrumental in our understand- 
ing of hadronic matter if proven correct. In particular, 
the HBT puzzle suggests we lack understanding of the 
last stages of the fireball evolution. Non-equilibrium is 
useful in this respect, since it affects both system volume 
0,|2^| and emission time [Ls| . 

The reason both equilibrium and non-equilibrium are 
compatible with data is that in a fit to yields the non- 
equilibrium phase space occupancies 7 S and "f q corre- 
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late with freeze-out temperature 21], making a distinc- 
tion between a T — 170 MeV equilibrated freeze-out 
(j q = 1,7 S < 1) scenario and a supercooled scenario 
where 7 9jS > 1 problematic. 

A related ambiguity is the difference between chemi- 
cal freeze-out (where particle abundances are fixed) and 
thermal freeze-out (where particles cease to interact). 
Equilibrium models generally assume a long phase be- 
tween these two points, which would alter considerably 
the multiplicity of directly detectable resonances. In 
a Non-equilibrium supercooled freeze-out, on the other 
hand, it is natural to assume that particle interaction af- 
ter emission is negligible. Once again, a reliable way to 
probe the extent of the reinteraction would be instrumen- 
tal for our understanding of how the fireball produced in 
heavy ion collisions breaks up. 

We have recently shown [14j that event-by-event fluc- 
tuations can be used to solve both of the dilemmas dis- 
cussed above. From the equations in the previous sec- 
tion, it is clear that the dependence of fluctuations of 
T and 7 g is different, allowing us to decouple these two 
variables. A higher temperature tends to decrease fluctu- 
ations w.r.t. the Poisson value expected from Boltzmann 
statistics (a = 1) , since it introduces greater correlations 
due to an increased contribution of resonances (D12 in 
Eq. |3J) . Vice-versa, increasing j q will rapidly increase 
fluctuations of quantities related to pions, due to the fact 
that at j q > 1 Xjr rapidly approaches e m ^l T , giving fluc- 
tuations an extra boost w.r.t. yields as per Eq. [§] 

In addition, comparing fluctuations to directly de- 
tected resonances probes the interval between chemi- 
cal and thermal freeze-out. Consider, for example, the 
K + /tt~ fluctuation. The top and the bottom terms in 
this ratio are linked by a large correlation term due to the 
if* (892) decay. This correlation probes the K*° abun- 
dance at chemical freeze-out, since subsequent rescat- 
tering/regeneration does not alter the fact that the K*° 
decay produced a n + and a K~ . On the other hand, a 
direct measurement of the K* / K~ ratio through invari- 
ant mass reconstruction measures the K* abundance at 
thermal freeze-out, after all rescattering. Hence, compar- 
ing the K/ir fluctuation to the K*/K~ ratio provides a 
gauge for effect of the hadronic reinteraction period on 
particle abundances. In particular, a sudden freeze-out 
model should describe both with the same parameters. 

The main experimental problem with fluctuation mea- 
surements is the vulnerability to effects resulting from 
limited detector acceptance. This difficulty can be less- 
ened, to some extent, by considering "dynamical" fluc- 
tuations, obtained by subtracting a "static" contribution 
which should be purely Poisson in an ideal detector. The 
exact nature of the subtraction varies from collaboration 
to collaboration, for example [27| using 

V dVn = ^ 2 -° 2 stat (11) 

(T x tn. t. , usually obtained through a Mixed event approach 
[28j, includes a baseline Poisson component, which for a 
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FIG. 2: Temperature and 7 9 statistical significance (Ptrue) 
profiles for the fit in Fig. 1. At each point in the abscissa a 
fit is performed varying all fit parameters except the one on 
the abscissa. Ptrue is calculated from the best fit parameters 

ratio N1/N2 can be modeled as 

«lat = J^) + J^ (12) 

as well as a contribution from detector efficiency and 
kinematic cuts. Provided certain assumptions for the 
detector response function hold (see appendix A of [2Sj]h 
subtracting a s t a t from a should yield a "robust" detector- 
independent observable. 

When using measured a dyn s in fits to experimental 
data, the data-sample should include a dyn , particle yields 
(which are needed to determine the Poisson contribution 
to a s t ci t as per Eq. I12f) and particle ratios. 

We have performed a fit incorporating all ratios from 
0, with the exception of the A++/p i _which the STAR 
collaboration has taken back, and the fi/fi, measured by 
STAR to be above unity (The latter measurement is a 
priori incompatible with statistical production \2^). In 
addition, we have included the preliminary value for the 
K/ir event- by-event fluctuation measured by STAR 27], 
as well as the published yield for cj) [3(J and n~ |3l| . 

Our fit parameters include the normalization (hope- 
fully related to the system "volume" at chemical freeze- 
out), temperature, X q , s ,i 3 and 7 9jS . We also require, by 
implementing them as data-points, strangeness, charge 
and baryon number conservation 
«s - s) = (Q) / (B) = «Q) / (B)) Au = OA.). 

The result is shown in Fig. ^ The fit gives an adequate 
description of all particle yields, including the resonance 
K* /K~ and A(1520)/A. It can also, adequately describe 
the event-by-event fluctuations of of K + /ir + or K~ /ir~ , 
but not the difference between them. Describing the dif- 
ference within the SHM is impossible without introduc- 
ing a sizable isospin fugacity A/3 ~ 0.96 which is in turn 
excluded by ratios such as ir + /ir~ and Q/B. We wait 
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FIG. 3: Sensitivity of the /T/7T fluctuation io^?^) to 7 9 and 
temperature for the fit in Fig. 1. At each point in the abscissa, 
a fit is performed (only yield data-points used) varying all fit 
parameters except the one on the abscissa, is calculated 

using the best fit parameters, dot-dashed lines refer to the 
experimental limits |27| for the a K + /„+ (blue) and a K -i^- 
(red). Dotted vertical lines give published non-equilibrium 
best fit parameters |2l|. while the long-dashed line and shaded 
area gives published equilibrium best fit |2ll 124 125||. 

for published results, and systematic error computation 
before studying this discrepancy further. 

Fig. |3 shows that the correlation between T and j q 
disappears when fluctuations are taken into account, and 
that the fit tightly constrains 7 9 well above the equilib- 
rium value, forcing temperature to be about 140 MeV, 
in good agreement of the prediction of the supercooled 
hadronization scenario. In accordance with this scenario, 



our model can describe both the K/ir fluctuation and the 
directly observed K* /K~ value with the same statistical 
parameters, as is expected for a system where hadronic 
reinteraction (the phase between chemical freeze-out and 
thermal freeze-out) is negligible. 

Fig. [^confirms that fluctuations arc indeed responsible 
for the tight constraining of T and j q . The best fit at light 
quark equilibrium misses cr^™ by several standard devi- 
ations. Introducing exact conservation for strangeness 
within the observed window would increase the discrep- 
ancy by reducing ok HU- It is only through 7 9 > 1 
that cr^™ increases until it becomes compatible with the 
experimental value. 

In conclusion, we have used preliminary experimental 
data to show that, at 200 GeV, the equilibrium model is 
unable to describe both yields and fluctuations within 
the same statistical parameters. The non-equilibrium 
model, in contrast, succeeds in describing almost all of 
the yields and fluctuations measured so far at RHIC, 
with the parameters expected from a scenario where non- 
equilibrium arises through a phase transition from a high 
entropy state, with super-cooling and oversaturation of 
phase space. We await more published data to determine 
weather the non-equilibrium model is really capable of 
accounting for both yields and fluctuations in all light 
and strange hadrons produced in heavy ion collisions. 

Work supported in part by grants from the U.S. De- 
partment of Energy (J.R. by DE-FG02-04ER41318) the 
Natural Sciences and Engineering research council of 
Canada, the Fonds Nature et Technologies of Quebec. 
G. T. thanks the Tomlinson foundation for support, and 
the QM2005 organizers for financial support during the 
conference. S.J. thanks RIKEN BNL Center and U.S. 
Department of Energy [DE-AC02-98CH10886] for pro- 
viding facilities essential for the completion of this work. 



[1] E. Fermi, Prog. Theor. Phys. 5, 570 (1950). 

[2] I. Pomeranchuk, Proc. USSR Academy of Sciences (in 

Russian) 43, 889 (1951). 
[3] L. D. Landau, Izv. Akad. Nauk Ser. Fiz. 17 (1953) 51. 
[4] R. Hagedorn, Suppl. Nuovo Cimento 2, 147 (1965). 
[5] S. Jeon, V. Koch, K. Redlich and X. N. Wang, Nucl. 

Phys. A 697, 546 (2002) 
[6] S. Jeon and V. Koch, Phys. Rev. Lett. 83, 5435 (1999) 
[7] S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000) 
[8] C. Pruneau, S. Gavin and S. Voloshin, Phys. Rev. C 66, 

044904 (2002) arXiv:nucl-ex/0204011 . 

[9] J. Adams et al. [STAR Collaboration], Phys. Rev. C 68, 

044905 (2003) 

[10] K. Adcox et al. [PHENIX Collaboration], Phys. Rev. 

Lett. 89, 082301 (2002) 
[11] V. V. Begun, M. Gazdzicki, M. I. Gorenstein and 

O. S. Zozulya, Phys. Rev. C 70, 034901 (2004). 
[12] G. Torrieri, et al., Comm. in Computer physics in press 
|arXiv:nucl-th /0404083 



[13] J. Cleymans, K. Redlich, Phys. Rev. C 60, 054908 (1999). 
[14] G. Torrieri, S. Jeon, J. Rafelski Submitted to PRL 

arXiv:nucl-th/0503026 . 
[15] J. Letessier and J. Rafelski, Cambridge Monogr. Part. 

Phys. Nucl. Phys. Cosmol. 18, 1 (2002), 
[16] See most advanced physical chemistry textbooks, or 

http:/ /www2. mcdaniel. edu/ Chemistry /ch3 07. notes/ 

Chemical%20Equilibrium.html 
[17] J. Rafelski and J. Letessier, Phys. Rev. Lett. 85, 4695 

(2000). 

[18] T. Csorgo and L. Csernai, arXiv:hep-ph/9312330 

[19] J. Letessier and J. Rafelski, Int. J. Mod. Phys. E 9, 107, 

(2000), and references therein. 
[20] J. Rafelski, J. Letessier, J. Phys. G 30, S I (2004). 
[21] J. Letessier and J. Rafelski, arXiv:nucl-th/0504028 
[22] D. Flierl et al. [NA49 Collaboration] , 

|arXiv:nucl-ex/041004l| 
[23] J. Rafelski, J. Letessier and G. Torrieri, Phys. Rev. C 72, 

024905 (2005) |arXiv:nucl-th/0412072] . 



5 



[24] O. Bar annikova [STAR Collaboration], 

arXiv:nucl-ex/0403014 and references therein. 
[25] P. Braun-Munzinger, D. Magestro, K. Redlich and 

J. Stachel, Phys. Lett. B 5 18, 41 (2001). 

[26] G. Torrieri and J. Rafelski, arXiv:hep-ph/0409160 
[27] S. Das [STAR Collaboration] , |arXiv:nucl-ex/0503023] 
[28] C. Pruneau, S. Gavin and S. Voloshin Phys. Rev. C 66, 
044904 (2002) arXiv:nucl-ex/0204011 . 



[29] M. Bleicher et al., Phy s. Rev. Lett. 88, 202501 (2002) 

arXiv:hep-ph/0111187 . 
[30] J. Adams et al. [STAR Collaboration], Phy. Lett. B in 

press, arXiv:nucl-ex/0406003 . 
[31] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 

92, 112301 (2004) |arXiv:nuc l-ex/0310004 . 



